Magnetic Field Effect on the Supercurrent of an SNS junction 
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In this paper we study the effect of a Zeeman field on the su- 
percurrent of a mesoscopic SNS junction. It is shown that the 
supercurrent suppression is due to a redistribution of current- 
carrying states in energy space. A dramatic consequence is 
that (part of the) the suppressed supercurrent can be recov- 
ered with a suitable non-equilibrium distribution of quasipar- 
ticles. 

PACS numbers: 74.50. +r, 74.80. Fp, 85.25.-j 

Recently, there has been a revival of interest in the 
proximity effect Q which occurs when normal metals (N) 
are placed in contact with superconductors (S). This is 
due to the availability of submicron fabrication technol- 
ogy and low temperatures. One of the main lessons learnt 
in the study of these hybrid structures is that, at suffi- 
ciently low temperatures, due to the suppression of inelas- 
tic scattering etc., it is essential to properly understand 
the contribution from individual energies, |2j |5[]. This ap- 
plies to, e.g., the conductance between N and S in contact 
as well as supercurrent in an SNS junction 

For the latter systems, it is in particular useful to un- 
derstand the supercurrent being carried at each individ- 
ual energy, i.e., the current-carrying density of states 
JVj(e). This quantity is the ordinary density of states 
weighted by the current that the states carry. The super- 
current is the integral over the energy e of the product of 
Nj(e) and the occupation number n(e). 

Here we consider again an SNS junction but with an ap- 
plied magnetic field. It is well-known that magnetic field 
in general suppresses the supercurrent. This can arise 
from two completely different mechanisms JTo| ] . First, it 
can be due to the coupling of the magnetic field to super- 
conductivity via the vector potential. I shall not discuss 
this effect here. This effect can be suppressed in geome- 
tries where the area perpendicular to the magnetic field 
is sufficiently small. The second mechanism is due to the 
Zeeman energy. Since pairing is singlet in s-wave super- 
conductors, a physical picture commonly used is that the 
Zeeman field is a pair-breaking perturbation and hence 
the Cooper pair amplitude decays in space faster in the 
presence of the field. In a dirty normal metal with diffu- 
sion coefficient D, this decay length ~ y/D/h [fioj for an 
energy splitting h. For an SNS junction, the supercurrent 
thus decreases with field due to a reduction of coupling 
between the two superconductors. 

This, however, is not entirely the full picture. As will 
be shown below, the main effect of the Zeeman splitting 
is to shift the current-carrying density of states in energy 



space. This is analogous to the behavior of the ordinary 
density of states under h [jl(J. The pairing correlation 
between the two superconductors remains long-ranged at 
appropriate energies. The supercurrent decreases ( c./., 
however, below) because of this mentioned shift and the 
associated change in the occupation of the states (see 
below for details). 

A dramatic consequence of the above is that, un- 
der a suitable non-equilibrium distribution of quasipar- 
ticles, one can recover the suppressed supercurrent. I 
shall demonstrate this using an experimental arrange- 
ment studied in ref A suitable applied voltage can 
recover (partly) the supercurrent suppressed by a mag- 
netic field. 

Consider thus a quasi-one dimensional dirty metal wire 
(N') of length L connecting two superconductors. We 
shall always assume that the junction is in the dirty limit. 
The supercurrent carrying density of states can be ob- 
tained via the angular averaged retarded green's function 
matrix g. For the present purposes it is sufficient to con- 
sider the retarded component and I shall leave out the 
usual superscript R. g obeys the normalization condition 
g 2 = — 7r 2 l and the Usadel equation. The latter reads, 
for position x within N' (0 < x < L) and a magnetic field 
B along the z direction, 



[ef 3 + ha 3 ,g] + —d x (gd x g) = 

7T 



(1) 



Here e is the energy with respect to the Fermi level, and 
h = [i e B. In order to avoid confusions I shall pretend 
that electrons have a positive magnetic moment and thus 
identify the directions of the magnetic moment and spin, 
with up (down) being the states with lower (higher) Zee- 
man energy, g at the boundaries x — and L are given by 
its corresponding values for the equilibrium superconduc- 
tors. I shall assume that the magnetic field is perfectly 
screened in S. In this case the boundary conditions at the 
superconductors are given by 



.9 = — 71 



ef 3 - A 



(2) 



with suitable gap matrices A reflecting the phase differ- 
ence x between the two superconductors. 

Eqn (|l|) can be simplified by noting that it is block- 
diagonal, since the pairing is singlet. In the usual 4x4 
notation, the elements associated with the 1st and 4th 
rows and columns are decoupled from those of 2nd and 
3rd (as already noted in, e.g., ref 11 12]) Moreover, the 
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matrix equations for these submatrices have the same 
structure as that in zero field except e — > e ± h, corre- 
sponding to magnetic moment parallel and antiparallel 
with the external field. It is then convenient to introduce 
separately the current-carrying density of states for each 
spin direction: 

NJ(e)=<p x N°(p,e,x)> (3) 

where N a (p, e, x) is the density of states for spin er ( = f 
or |), p the direction of momentum. Here I have chosen 
to label the states with the spin direction of the parti- 
cles. Note that, e.g., the up-spin particles are associated 
with the down-spin holes (c./. above). For a given spin 
a, Nj(e) is related to the appropriate sub-matrix of the 
green's function g via formulas analogous to those in zero 
field (see, e.g. §) 

The total (number) current density at T = is given 
by the integration of over the occupied (negative) 
energy states: 

Js=v f I de[Nl(e)+N l j(e)} (4) 

J — -OO 

We shall also introduce NJ v (e) = |[iV)(e) + iVj(e)] as 
the spin averaged current-carrying density of states. J s 
is related to NJ"(e) by 

J s = 2v f [ de[NT(e)}, (5) 

J — OO 

exactly the same formula as in zero field. Here Vf is the 
Fermi velocity. 

We shall confine ourselves to the case of long junctions 
( Ed « |A| ), and for defmiteness choose |A| = IOOEd- 
Here Ed = D/L 2 is the Thouless energy associated with 
N'. The behavior of the current-carrying density of states 
Nj(e) in zero field has already been studied in detail in 
HH . I shall only mention some of the more relevant fea- 
tures below. Nj(e) vanishes for all energies at % = 0. A 
typical case for other phase differences < x < f ( we 
shall always restrict ourselves to this range, the other 
cases can be obtained by symmetries) is as shown by full 
line in Fig [l| Nj is odd in the energy variable e. Its 
major feature consists of a positive peak (labelled by + 
in Fig P at energies of several times the Thouless energy 
Ed below the fermi energy, and a corresponding negative 
peak above (labelled by — in Fig. [j]). Also seen are the 
small undulations as a function of energy for larger ener- 
gies. This oscillatory behavior is a result of the difference 
in wave-vectors for the participating particles and holes 
undergoing Andreev reflection at given e. For one dimen- 
sion and in the clean case the pairing amplitude / (the 
off-diagonal elements of g in particle-hole space) oscillates 
as e ± 2l ( e / v f) x . i n the present dirty three dimensional case 
the same physics results in e.g., f ~ expi [(1 — i)yJt/Dx] 
for the linearized Usadel equation (i.e. the limit of small 



pairing amplitudes) At a fixed position, the pairing am- 
plitude and hence the coupling between the two super- 
conductors oscillates as a function of energy. Nj{e) also 
vanishes for |e| below a few times Ed, where the ordinary 
density of states also vanishes. Both the magnitude of 
this 'minigap' and the position of peaks mentioned above 
decrease with increasing phase difference, vanishing as 
X — > tt (where Nj itself also vanishes for all e's). 

The behavior of Nj(e) under a finite field is also 
shown in Fig |f], where I have chosen an intermediate h 
( Ed « h « |A|) for clarity. As mentioned it is con- 
venient to discuss the current-carrying density of states 
separately for each spin direction under the presence of 
h. For magnetic moment along the applied field, the cur- 
rent carrying density of states is roughly (c.f. below) 
that of zero field except shifted in energy by — h. i.e. 
Nj(e) » AjT°(6 + h). Correspondingly N l j for magnetic 
moment pointing in the opposite direction is shifted up 
in energy. At fields h not too small compared with A, 
there is a correction to this picture because, if we assume 
perfect screening of the magnetic field inside the super- 
conductor as we are doing, the replacement e — ► e ± h 
in the Usadel equation does not apply for the boundary 
condition (||) at x = and L. This correction is negligible 
if h << |A| and increases with increasing h. 

In this picture, the reason that the supercurrent at fi- 
nite h is suppressed (in general) from that of zero field 
is not pair-breaking, at least for h << |A|. Rather, it 
is because at finite field ~ Ed, some of the states that 
have positive contributions to J s were orginally occupied 
at h = but are now empty (+ J.)., whereas some which 
were orginally empty are now occupied ( — T) an d con- 
tribute a negative current. 

By the above reasoning, the presence of the magnetic 
field has a non-trivial effect on the current-phase relation- 
ship. An example for a small h << |A| is shown in Fig 
^. In zero field I(x) is roughly like a sine function except 
for a small tilt towards x ~ n - When h increases from 
zero, one sees that I s first starts to decrease for x near 
7r while I s at smaller x is unaffected. Only at larger h 
would I s begin to be suppressed there. This can be read- 
ily understood by considering the behavior of Nj under h 
discussed before. Recall that at zero field Nj has a mini- 
gap of order Ed but x dependent, being smallest when x 
is near w. Thus when the field h is increased from 0, I s 
at larger x would be suppressed first since at these x' s j 
a smaller h is needed to shift the antipeak — f (the peak 
+ i) to below (above) the fermi level. 

At higher h and near \ ~ n t he current also oscillates 
with x- (See, in particular h = 6 in Fig. [|, where I s 
becomes negative for x slightly less than n, vanishing 
again at x = 7r )- These features are due to the undulatory 
structure of N] as a function of e (see Fig |l|). At these 
higher fields the weaker bumps and troughs of Nj (not 
labelled) cross the fermi level successively. They do so 
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at fields which are \ dependent. Their amplitudes also 
depends on \- 

Since the major effect of the magnetic field is not a 
suppression of Nj (except for /i^|A|) but a redistribu- 
tion in energy space, the supercurrent suppression by the 
magnetic field can, to a certain extent, be recovered by a 
suitable distribution of quasiparticles. Here we consider 
the case of a 'controllable Josephson Junction', studied 
in experimentally and ||[| theoretically. The device 
configuration is shown schematically in the inset of Fig 
||. The two superconductors S, in general with a phase 
difference x, are at zero voltage. Equal but opposite volt- 
ages V are applied on the normal (N) reservoirs. The N 
and S reservoirs are connected by quasi-one dimensional 
normal wires as shown. We are interested in the effect of 
the voltage V on the supercurrent I s between the S reser- 
voirs. At V = 0, the distribution function n(e) is of the 
usual equilibrium form and is given by 1 for e < and 
otherwise (T = 0), as shown by dotted lines in Fig |[ Un- 
der a finite V, n(e) for |e| < eV becomes 1/2 [Q (full line 
in Fig ||), i.e. its effect is to transfer half of the quasipar- 
ticles for — eV < e < to the region < e < eV. In zero 
field such a non-equilibrium distribution in general leads 
to a decrease of the supercurrent between the S reservoirs 
(see [^-||), since usually this corresponds to decreasing 
(increasing) the occupation of states which contribute a 
positive (negative) current [ + (— ) in Fig [jj (full line in 
Fig U p5| ; see ref for the oscillations at higher V's) 

The effect of V on I s for finite h is also shown in Fig |]. 
At the fields chosen (>> Ed) I s at zero voltages have es- 
sentially decreased to zero (see also below). As claimed, 
for a given h, a suitable choice of V may 'enhance' the su- 
percurrent. For the present case where Ed « h « |A|, 
this enhancement is particularly spectacular for eV ~ h. 
These features can be understood by examining the spin- 
averaged current-carrying density of states Nj v (e), also 
shown in Fig |3|. Nj v (e) is odd in energy. The behavior 
of Nf follows directly from N] 1 in Fig |[ For a given 
h 0, when V is increased from zero, we have a trans- 
fer of the particles as described before from the states 
near the antipeak just below the fermi level ( due to — f) 
to the peak above (arising from + |). The supercurrent 
thus increases. The strongest enhancement of the super- 
current occurs at eV ~ h, since the region where this 
transfer occurs just covers the antipeak below the fermi 
level. For Ed « h « |A|, roughly half of the super- 
current at zero voltage can be recovered at eV ~ h. This 
is because at this voltage, the contributions from the re- 
gion — eV < e < eV cancel among themselves and we 
are left with the integral over states with e < — eV. The 
integral of Nf over the energy in this region is roughly 
half that of zero field, (compare Fig || and Fig [j]) 

The above statements are not quantitatively precise 
due to: (i) Nj are not simple shifts of Nj in energy and 
(ii) oscillatory structures of , both mentioned before. 



For V — 0, as can be seen in Fig || the supercur- 
rent does not simply decay monotonicially with increas- 
ing h but rather shows a damped oscillation. Such a 
behavior has already been pointed out before |lT| . |l2| . 
and was explained in terms of the oscillation of the pair- 
amplitude with h at a given distance (the separation be- 
tween the two superconductors). The present work pro- 
vides a slightly different but closely connected perspec- 
tive. I s oscillates with h because Nj does so as a function 
of energy. The field h shifts NJ in energy space. The 
undulatory behavior results when regions of alternating 
signs of Nj shift through the fermi level. 

Some of the physics discussed in this paper, such as the 
effect of h on the current-phase relationship, is applicable 
beyond the dirty limit. I shall however defer these to a 
future study. 

There is recently strong interest in the physics of su- 
perconductors in contact with a ferromagnetic material. 
Many papers have simply modelled the ferromagnet F 
with a Stoner field, e.g. jl6| (c./., however |l7[]) Within 
this model the Stoner field is formally equivalent to the 
Zeeman field here. A much discussed topic is the effect 
of this field h on the number of conduction channels jl6| . 
This effect is important only when h is comparable to 
the fermi energy of N' and has not been included in the 
present calculations. 

Part of the paper was written using the facilities of the 
Chinese University of Hong Kong. I thank Profs. P. M. 
Hui and H. Q. Lin for their help. 
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FIG. 1. Nj (dot-dashed and dashed respectively and 
in units of N f l/6L Q) for a junction with A = W0E D 
at h = 3O-E15. Also shown is Nj(e) in zero field (full line). 




FIG. 2. Current-phase relationships for magnetic 
fields h displayed in the legend, h is in units of Ed, Is in 
units of Ed/2Rm- Here Rm is the normal state resistance 
of N' = 2Nf DSL^ 1 with S the cross-section area of 
N'). 
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FIG. 3. Explanation of the effect of a voltage 
on the supercurrent for the device shown in the inset. 
The occupation numbers n(e) (left scale) are shown for 
V = (dotted) and V = 30 (full line). Also shown is 
the spin-averaged current-carrying density of states Nj v 
(right scale), e and eV are in units of Ed- 
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FIG. 4. Supercurrent as a function of voltage V be- 
tween the two superconductors S of the device depicted 
in the inset of Fig |!| The values of h, in units of Ed, are 
shown in the legend, x = t/2. 
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